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We apply our method of complex scaling, valid for a general class of potentials, in a search for 
nucleon-nucleon S-matrix poles up to 2 GeV laboratory kinetic energy. We find that the realistic 
potentials JISP16, constructed from inverse scattering, and chiral field theory potentials N 3 LO and 
N 2 LO op t support resonances in energy regions well above their fit regions. In some cases these 
resonances have widths that are narrow when compared with the real part of the S-matrix pole. 
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Introduction. There is a long-standing interest and 
considerable recent progress in the theoretical charac¬ 
terization of nuclear resonant states. A resonant state 
is fully characterized by its position in the energy plane 
and its width, which determines how fast the state will 
decay. One could in general solve the time-dependent 
Schrodinger equation to study the characteristics of res¬ 
onant states [lH3(, which is a demanding computational 
process. On the other hand, the time-independent many- 
body methods that deal with the description of resonant 
states in nuclei are under development and exhibit ap¬ 
pealing computational features. These time-independent 
methods can be divided into real energy and complex 
energy approaches. 

The spectrum of a real nuclear Hamiltonian consists 
of negative and positive energy states. While the nega¬ 
tive energy spectrum is discrete (bound states), the pos¬ 
itive energy spectrum may have a richer structure with 
resonant states among scattering or continuum states. 
Hence real-energy approaches require criteria for identi¬ 
fying a resonant structure and for assigning a position 
and a width to them. 


In the domain of C 2 integrable basis expansion meth¬ 
ods this is usually achieved through C? stabilization 
methods Si or methods that evaluate the real Con¬ 
tinuum Level Density (CLD) [§]. The CLD usually pro¬ 
duces an approximate Breit-Wigner distribution in the 
region of the resonant state whose parameters could be 
determined by a fit. The CLD method has been success¬ 
fully applied to atomic systems (zj, nuclear clusters Q 
and in mean-field approaches for describing quasiparticle 
resonant states [§]. 


The name stabilization, arises from the fact that one 
does not need the knowledge of the asymptotic wavefunc- 
tion in order to determine the resonant parameters. On 
the reaction side, based on R-matrix considerations Ha¬ 


ul , and assuming the single channel approximation, res¬ 


onant parameters can be determined by the behavior of 
the phase-shift function of energy 5(E) around the reso¬ 
nant position; in particular the position is defined as the 
inflection point of 5(E) (maximum energy derivative of 
5(E)) and the width is defined as T = dS J dE \ E=E r , where 


E r is the inflection point. Such formulas where employed 
recently in microscopic R-matrix calculations [13] to ex¬ 
tract widths from realistic nucleon-nucleus phase-shifts. 
Though the R-matrix parametrizations have been very 
successful, the formulas become less transparent in the 
multi-channel case and when they are applied for the 
description of broad resonances (see discussion in HI)- 
Furthermore, for broad resonances, R-matrix analyses 
become more dependent on channel radii and boundary 
conditions [3]. Finally, combining formulas and assump¬ 
tions from different theories/models for the calculation of 
an observable increases the possibility of uncontrollable 
errors. 

The complex energy formalism serves as a potentially 
fruitful alternative for the characterization of the reso¬ 
nant parameters. It has been shown that once the R- 
matrix, S-matrix and T-matrix are analytically contin¬ 
ued to the complex energy plane, the extraction of reso¬ 
nant parameters becomes independent of boundaries and 
radii [is]. Apart from this practical issue, some physi¬ 
cal phenomena may have a more natural interpretation 
once the theory is developed in the complex energy plane 
(e.g. thermo-nuclear reactions fl7l fl8|). In the complex 
energy formalism, the Gamow (resonant) states, i.e. the 
solutions of the Schrodinger equation which satisfy purely 
outgoing boundary conditions (complex wave number k), 
play a dominant role. It was shown by Berggren [19] that 
resonant states, when accompanied by non-resonant con¬ 
tinuum states, form a complete set, an important prop¬ 
erty that gives rise to Berggren basis expansion meth¬ 
ods either in a Configuration Interaction (Cl ) framework 
[ifll - Eiij . Coupled Cluster framework [ifil - [ 28 l ] or reaction 
theory framework [29l435| . Expressing the Hamiltonian 
in such a complex energy, orthonormal non C 2 integrable 
basis, automatically allows its spectrum to support res¬ 
onant and also non-resonant continuum states. In addi¬ 
tion, when the Berggren basis is used in a reaction frame¬ 
work the detailed knowledge of the boundary condition 
at large distances is not crucial. 

The Complex Scaling (CS) method also belongs in the 
category of complex energy formalisms. The Aguilar- 
Balslev-Combes (ABC) theorem 36, 37] establishes that 
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once the Hamiltonian coordinates are rotated, the reso¬ 
nant states are independent of the rotation and behave 
asymptotically as bound states. Consequently, one could 
use the technology that has been established for bound 
states in order to describe resonant and scattering phe¬ 
nomena. Furthermore, the CS method has been success¬ 
fully applied in nuclear physics 38 -42| (see also [ 4 !] for 
an application of CS in a deformed nuclear mean-field). 
We recently showed [Hj that this method may be applied 
to the most general cases of non-local nuclear potentials. 

In this work we apply the CS method to nucleon- 
nucleon (NN) scattering spanning the range from thresh¬ 
old to 2 GeV laboratory kinetic energy, which exceeds the 
fitting range of most NN potentials. We elect to retain 
non-relativistic kinematics throughout as the interactions 
are derived for a non-relativistic scattering framework. 
We employ three different realistic NN interactions and 
we find resonant poles at laboratory kinetic energy of 
about 600 MeV, or at about 2.2 GeV in the total cen¬ 
ter of mass energy. Some of these poles correspond to 
narrow resonant states. 

According to the SAID data analysis group 45-471] 
(see also [48() there exist resonance-like structures, poles 
of the S-matrix, in the 1 D 2 , 3 F 3 uncoupled and coupled 
3 P 2 - 3 F 2 channels. Recently a resonant-like structure was 
also found by the WASA-at-COSY collaboration and the 
SAID analysis group in the 3 D 3 - 3 G 3 coupled channel 
©ii Our CS calculations, in addition to showing 
resonances in these channels, also reveal resonance-like 
structures in the 3 Po and 3 Pi channels. We searched 
other channels up to and including L=4 without any ad¬ 
ditional signals of resonance-like structures. 

The study of dibaryon resonances could shed light on 
the reaction mechanism and aid in the interpretation of 
excited nucleonic states. It is also valuable to pin down 
the properties of dibaryon resonances as a potential link 
between Quantum Chromodynamics, hadron models and 
traditional low energy nuclear physics. In the work of 
Fibl 1tt| the resonant-like structures where identified by 
analytically continuing the T-matrix of the available data 
in the complex energy plane. Our goal is to simply iden¬ 
tify resonant structures with the CS method but not to 
study in depth the characteristics of the NN scattering at 
intermediate energies, something that would require the 
use of NN interactions that fit scattering data at higher 
energies, such as CDBonn 5l| or AV18 [52 L Such in- 
depth studies would be done relativisticly 53| and by 
treating properly A and/or Roper resonances degrees of 
freedom (see for example 0). Furthermore, we will 
not provide information on the possible decay paths that 
the resonant structures will follow, since we do not con¬ 
sider couplings to inelastic channels such as, NN —» 7rd 
or NN-> AN etc. In addition, the interactions we use, 
are modeling the short-range (high-energies) NN sector 
in different approaches and are fitted at lower labora¬ 
tory energies (< 350 MeV). Hence, we are not aiming 


at making predictions for the existence or absence of 
broad dibaryonic states. Nevertheless, it is worthwhile 
to discover that the NN interactions we employ, sup¬ 
port high energy resonant-like states above the A pro¬ 
duction threshold (1232 MeV). The existence of these 
resonances indicates such degrees of freedom suggest it 
may be important to examine what effect the explicit 
inclusion of the related degrees of freedom will have on 
finite nuclei. Simply stated, our goal at this point is 
to demonstrate that the CS method locates these reso¬ 
nances using three different realistic NN interactions in 
the conventional non-relativistic framework. 

Method and results. We apply the CS transformation 
to our Hamiltonian which consists of the relative kinetic 
energy T and the realistic NN interaction V between 
the nucleons. The complex rotated Hamiltonian has the 
form: 

H(r, 9) = e~ 2ie T + V(re ie ), (1) 


where 9 is the real CS rotation parameter. Nuclear po¬ 
tentials are usually non-local in coordinate space but, for 
simplicity and compactness, the potential is taken to be 
local in EmO]). In 44] we presented our methodology 
(see also [55[) for applying the CS transformation to a 
non-local potential which we also follow here. The time- 
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FIG. 1. (Color online) Spectrum of the CS Hamiltonian JT]) 
for the N 3 LO (a), N 2 LO op t (b) and JISP16 realistic potentials 
for the 3 Pi channel (or 1“ state) of the np system. There is a 
state that is invariant of the rotation angle 9 (rad), indicated 
by an arrow. In panel (d) we show a magnification of the 
N 3 LO Hamiltonian spectrum in the region of the resonant¬ 
like state at a sequence of cutoffs in the number of radial 
nodes (N) in the basis. This is the ^-trajectory. Energies are 
in the Center of Mass (CoM) frame. 

independent non-relativistic Schrodinger equation then 
becomes: 


H(r , 0)tf(r, 9) = E{6)^(r, 9). (2) 

where E is the energy in the Center of Mass (CoM) frame 
here and throughout this work. To be more precise, 
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the rotated non-Hermitian Hamiltonian operator and the 
Hermitian one, are related through the formula: 

H(r,e) = U(0)H(r)U(9)-\ (3) 


where U{9) stands for the non-Unitary CS transforma¬ 
tion operator. It is then expected that any quantum me¬ 
chanical operator O will be transformed under ([3]) in the 


case of CS (e.g. see [56| for an application on the dipole 


operator). In order to calculate its expectation value, 
one could either use the transformed operator and calcu¬ 
late its expectation value between CS solutions (9 ^ 0.0), 
or evaluate the bare operator and calculate its expecta¬ 
tion value using the coefficients from the back-rotated, 
^-independent solution. In [42j it was shown that both 
ways are equivalent and, in addition, the back-rotation 
transformation was stabilized via a smoothing process. 

In order to solve Eq.© we assume that the solution is 
a linear combination of orthonormal Harmonic Oscilla¬ 
tor (HO) basis states and we solve a complex symmetric 
Hamiltonian eigenvalue problem by diagonalization. The 
spectrum of the Hamiltonian contains resonant (bound 
states, resonances) and non-resonant continuum states. 
According to the ABC theorem, once the resonant state is 
revealed it remains invariant under CS rotations, whereas 
the non-resonant continuum states follow an approximate 
29 path in the complex energy plane. This is the complex 
stabilization criterion that is used in CS for the identifi¬ 
cation of the resonant state. In practice, due to the trun¬ 
cation of the underlying basis, there is a small variation 
of the resonant position with 9. It is then a consequence 
of the complex virial theorem [HtJ that the resonant state 
will be the one that corresponds to the minimal change 
of the real part of the energy with respect to 9. The 
method is also known as 0-trajectory method and it is 
a common practice in CS applications (see for example 
II El El). We also apply this stabilization technique 


and we check convergence of our results as a function of 
the basis dimension and the variations with the rotation 
angle. 

In Fig[l]we present the spectrum of the complex scaled 
Hamiltonian for the 3 Pi (1~) channel of the neutron- 
proton (np) system, for the JISP16 [60] and the two 
chiral effective field theory interactions N 3 LO [61] and 
N 2 LO op t [Hi]. The HO basis was characterized by hu> = 
40 MeV for JISP16 and N 2 LO op t and by fiui = 28 MeV 
for N 3 LO. For the N 3 LO potential we varied the rotation 
from 0 = 0.1 rad to 0.2 rad, for N 2 LO op t from 0 = 0.14 
to 0.24 rad and for JISP16 from 0 = 0.2 to 0.3 rad. The 
step in the 0 discretization was 0.004 rad. As we start 
rotating the coordinates and momenta of the Hamilto¬ 
nian, solutions that initially inhabit the real-axis (0 = 
0.0) start moving inside the complex energy plane. Dur¬ 
ing the rotation, when a state crosses the ordinates of a 
pole of the S-matrix, it remains there and does not follow 
the rotation of the other non-resonant continuum states. 
In Fig|T]we see clearly that all interactions support one 


state which is almost invariant with respect to the CS 
rotation parameter. For the calculations in panels (a), 
(b) and (c) we used a HO basis that consisted of a maxi¬ 
mum of N=30 radial nodes. The resonant parameters of 
this state where identified by applying the 0-trajectory 
method. Any Hamiltonian after the CS transformation is 
applied, becomes non-Hermitian and one needs to locate 
the resonant state or the stationary point. The stationary 
point is the one for which the difference in energy with 
respect to the 0 variation is minimal. The 0-trajectory is 
shown in panel (d) of FigJTJ We see that for each basis 
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FIG. 2. (Color online) Same as Fig[T]but for the 3 Po channel 
(CP state) and for a HO basis size of N=30. 

size (denoted by N) the 0-trajectory is different. Never¬ 
theless, the variations become smaller by increasing the 
basis size and the results for N=26 and N=30 start to co¬ 
incide for 0 > 0.17 rad. At this point we also notice that 
results are practically converged at N=18 and differences 
appear only in the fifth and sixth significant digit. 

In addition to the 0-trajectory method, the b- 
trajectory method is employed EH] in CS calculations. 
This b-trajectory method uses the dependence of the re¬ 
sults on the length parameter of the basis, in our case the 
HO length or, equivalently, hjj. We do not analyze the 
b-trajectory method here. Here, we simply note that the 
deuteron ground state (g.s.) ( 3 Si- 3 Di coupled channel) 
is very well converged (independent of Hlo) at N = 30. 
We will perform detailed b-trajectory analyses in future 
applications of CS to resonances in finite nuclei. 

We gather a sample of our results for the 3 Pi state in 
Table |H The sensitivity analysis regarding the 0 depen¬ 
dence and the basis convergence was also performed for 
the other partial waves. 

In Fig[2]we present the CS spectrum for the 3 Po chan¬ 
nel. For a basis size of N=30 the converged resonant¬ 
like states are (350.859 -i78.689) MeV, (305.859 -i44.827) 
MeV and (251.252 -i84.454) MeV for the N 3 LO, N 2 LO opt 
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TABLE I. 3 Pi resonant parameters in MeV as a function of 
the HO basis size. The choice of the energy was decided by the 
0-trajectory method. The energies are in the CoM of the np 
system above the np threshold (E=0). To obtain the width, 
the formula T = -2Im(E) could be applied. For the np system 
the nucleon mass is taken to be 938.2 MeV/c 2 . 


N 

N 3 ] 

LO 

N 2 LO op t 

JISP16 

14 

(214.767 

420.048) 

(292.114 

449.609) 

(314.471 

484.250) 

18 

(213.758 

420.139) 

(292.358 

450.039) 

(314.475 

484.244) 

22 

(213.784 

420.131) 

(292.323 

450.023) 

(314.473 

484.244) 

26 

(213.781 

420.134) 

(292.327 

450.020) 

(314.473 

484.244) 

30 

(213.781 

420.134) 

(292.327 

450.020) 

(314.473 

484.244) 


and JISP16 interactions respectively. 

We now focus our attention on the N 3 LO and JISP16 
interactions, since they are both fitted to higher ener- 
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FIG. 3. (Color online) Spectrum of the 3 D 2 channel with the 
N 3 LO (a) and JISP16 (b) interactions. Notice the two poles 
in the case of JISP16. In panel (c) it is the coupled 3 D 3 - 3 F 3 
spectrum for the N 3 LO interaction. 


gies than N 2 LO opt . The 3 D 2 channel shown in Fig. [3] 
presents an interesting case since both interactions sup¬ 
port resonant states and, in the case of JISP16, we find 
two resonances. According to our calculations the N 3 LO 
resonant position is at (371.111 -i93.169) MeV. We no¬ 
tice that at an energy of about 400 MeV for the real part 
and about -150 MeV for the imaginary part, there is also 
a state that shows a stabilization pattern with respect 
to the rotation parameter. This state, however, is not 
as stable as the other cases we examine in this work, so 
we do not investigate it further. For the JISP16 interac¬ 
tion we observe clearly two resonant positions in the 3 D 2 
channel : (153.155 -i92.725) MeV and a broader struc¬ 
ture at (311.008 -i210.065) MeV (see panels (a) and (b) 
of Fig. [3]). 

In panel (c) of Fig[3]we show the spectrum for the cou¬ 


pled 3 D 3 - 3 G 3 channel for the N 3 LO potential which sup¬ 
ports a resonant state at a position: (219.218 -i75.162) 
MeV. For the same channel, not shown here, the JISP16 
interaction supports a broader structure at (200.084 - 
il 15.138) MeV. In this JISP16 case the rotation angle 
was varied from 9 = 0.32 to 0.42 rad in order to reveal 
the broader resonant position. 

In Fig. [4] we show our results for the uncoupled 3 F 3 
and the coupled 3 P 2 - 3 F 2 channels using the N 3 LO and 
JISP16 forces. For these partial waves we find multi¬ 
ple resonances supported by the interactions. The chiral 
potential supports at least four resonant-like structures 
for the 3 F 3 channel and for the coupled 3 P 2 - 3 F 2 chan¬ 
nel, where in the latter case we also see some indications 
for resonant structures at energies larger than 400 MeV 
CoM. According to our calculations, for the JISP16 in¬ 
teraction we find two resonant-like structures for each of 
the 3 F 3 and 3 P 2 - 3 F 2 channels. After applying the 9- 



FIG. 4. (Color online) Same as in Fig. [3] but for the 3 F 3 
channel (panels (a),(b)) and 3 P 2 - 3 F 2 (panels (c),(d)) for the 
N 3 LO and JISP16 interactions respectively. See also Table ITU 


TABLE II. Resonant parameters that correspond to the ar¬ 
rows of Fig. U computed with the 0-trajectory method for 
the uncoupled 3 F 3 and coupled 3 P 2 - 3 F 2 channels, using the 
N 3 LO and JISP16 interactions. 



N 3 

LO 



JISP16 


3 

f 3 

3 p 2 

- 3 F 2 

3 

f 3 3 p 2 - 

- 3 f 2 

(227.0 

4124.3) 

(247.8 

499.8) 

(242.3 

4173.4) (161.2 

495.2) 

(307.2 

4182.5) 

(318.9 

4120.5) 

(342.6 

4315.6) (342.4 ■ 

4213.2) 

(360.5 

4227.9) 

(353.5 

4160.6) 




(435.3 

4273.2) 

(392.0 

4179.1) 





trajectory stabilization method we gather the results in 
Table nn We notice that the states are broader than the 
structures we find for the 3 Po and 3 Pi channels. The 
positions of the resonant states in each channel depend 
on the interaction we use. In this sense, the situation is 
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distinctly different from the deuteron ground state that 
all the realistic NN interactions are constrained to de¬ 
scribe with high accuracy. Of course, this dependence 
of the scattering resonances on the NN interaction may 
be expected since we are investigating a kinematic region 
which is sensitive to short range nuclear physics, a sector 
that each NN interaction models differently. 

To show that the resonance locations are fixed by NN 
model assumptions and not by the couplings of low and 
high momenta alone, we made a test which demonstrates 
that the 3 Pi resonant state for the N 3 LO potential is 
invariant under similarity renormalization group (SRG) 
transformations [63|, down to a scale of A = 1.5fm _1 . A 
scale of A = 1.5fm , taken as a maximum relative mo¬ 

mentum, corresponds roughly to an energy of 187 MeV 
in the laboratory frame. The real part of energy of the 
3 Pi resonant-like state (see Table |U) is above the energy 
implied by this SRG scale (i.e. E l r a e b s ~ 2-214 = 428 MeV). 
Intuitively, one would think that the NN physics up to 
this SRG scale will remain invariant but above this scale 
one might expect some changes so we are motivated to 
test for SRG scale invariance using this resonance as the 
test case. 

Using the N 3 LO interaction in the 3 Pi channel we per¬ 
formed three SRG transformations, that reduce the cou¬ 
plings of low momentum to high momentum states and 
change the off-shell components of the potential, leaving 
the on-shell properties invariant. In Fig|5]we observe that 
the position of the 3 Pi resonant-like state is invariant 
under these SRG evolutions, for the SRG transformation 
scale as low as A = 1.5fm _1 starting from a bare N 3 LO 
potential. For these calculations the rotation angle was 
rotated from 9 = 0.1 rad to 0.2 rad with a step of 0.05 
rad. 

Conclusions. In this work we applied the CS method 
to study resonant features of the NN interaction in en¬ 
ergy regions above the usual 300 to 350 MeV laboratory 
energy. While this region exceeds the energies for which 
these non-relativistic interactions are developed, we ob¬ 
tain resonances motivating further exploration and we 
demonstrate the robust characteristics of our techniques. 

We analyzed the resonant features that we found in 
several channels and we studied the stability of the results 
using the 0-trajectory method whose validity has been 
demonstrated for local and schematic potentials. In our 
work we find similar convergence patterns for resonant 
states to what other authors have found. Among the 
numerous NN resonances that we found, the one in the 
3 Pi channel appeared close to the real energy axis. 

In our formalism, the widths of the states should be 
viewed as total decay widths, and since we did not con¬ 
sider special decay channels in our analysis, we cannot 
say if the np system will decay by either emitting mesons 
or baryon resonances. Our numerical results show that 
the positions of these states in the complex plane de¬ 
pend strongly on the form of the underlying interac- 



Re(E) (MeV) 

FIG. 5. (Color online) 3 Pi channel CS spectrum for the bare 
N 3 LO force and for three different SRG evolution scales. The 
resonant-like structure remains invariant. The green circle 
that corresponds to SRG A = 2.5 fm -1 practically overlaps 
with the results of the other SRG evolution scales and it is 
not visible. A solution is observed that corresponds to a A 
= 1.5 fm -1 with an imaginary part of about -250 MeV. This 
is a non-resonant continuum state which happened to depart 
from the approximate 29 line. 

tion. Hence these NN-dependent S-matrix poles are dis¬ 
tinguished from the deuteron pole whose properties are 
shared by all realistic NN interactions to high accuracy. 

We focused on a specific force in one channel, in partic¬ 
ular the chiral N 3 LO in the 3 Pi channel, to demonstrate 
that the position of the resonant-like state is invariant 
under the action of SRG transformations. 

Finally, it worth mentioning that analysis of NN scat¬ 
tering data has not, to our knowledge, reported resonant 
poles for the 3 Po and 3 Pi partial waves, but rather for 
the 1 D 2 , 3 F 3 , 3 P 2 - 3 F 2 and recently for the 3 D 3 - 3 G 3 . 

Among our achievements, we have demonstrated the 
stability of the CS method for identifying resonant states. 
This demonstration helps support the adoption of these 
techniques for investigating resonances in finite nuclei 
where we anticipate the interest will be in low-energy 
applications. 
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